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ABSTRACT 

Signs of damping wing absorption attenuating the Lyman-a emission line of the 
first known z ^ 7 quasar, ULAS J1120+0641, recently provided exciting evidence of 
a significantly neutral IGM. This long-awaited signature of reionization was inferred, 
in part, from a deficit of flux in the quasar’s Lyman-ct emission line based on predic¬ 
tions from a composite of lower-redshift quasars. The composite sample was chosen 
based on its C IV emission line properties; however, as the original study by Mort- 
lock et al. noted, the composite contained a slight velocity offset in C IV compared to 
ULAS J1120-I-0641. Here we test whether this offset may be related to the predicted 
strength of the Lyman-a emission line. We conhrm the significant 10 per cent at 
r.m.s.) scatter in Lyman-a flux for quasars of a given C IV velocity and equivalent 
width found by Mortlock et al. We further find that among lower-redshift objects 
chosen to more closely match the C IV properties of ULAS J1120-1-0641, its Lyman-a 
emission falls within the observed distribution of fluxes. Among lower-redshift quasars 
chosen to more closely match in C IV velocity and equivalent width, we find that ULAS 
J1120-I-0641 falls within the observed distribution of Lyman-a emission line strengths. 
This suggests that damping wing absorption may not be present, potentially weaken¬ 
ing the case for neutral gas around this object. Larger samples of z>7 quasars may 
therefore be needed to establish a clearer picture of the IGM neutral fraction at these 
redshifts. 

Key words: quasars: individual: J1120-I-0641 quasars: emission lines intergalactic 
medium cosmology: observations dark ages, reionization, first stars 


1 INTRODUCTION 


The reionization of hydrogen in the intergalactic medium 
(IGM) is believed to coincide with the buildup of the first 
galaxies. Determining when and how the IGM became ion¬ 
ized can therefore deliver unique insights into the earliest 
epochs of galaxy formation. Evidence from the Lyman-a 
forest in high-redshift quasars indicates that the process 
largely finished by z ~ 6 (Fan et al. 2006 Becker, Rauch 


|fc Sargent) |2007[) . While the exact timing remains unknown 


the recent Planck results (Planck Collaboration XIII 20151 


favour models of late reionisation, which motivates the ob¬ 
servational search for signs of reionisation. 

The recent discovery of the most distant quasar known. 


ULAS J1120-I-0641 (Mortlock et al. 20111, with a redshift of 


7.0842 ± 0.0004 (Venemans et al. 20121, provides a unique 


opportunity to probe the epoch of reionisation at a time 
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when the process was potentially just ending. As noted by 


Mortlock et al. (20111, this object exhibits strong absorp¬ 


tion on the blue side of its Lyman-a emission line. The red 
side of its Lyman-a line may also be somewhat weaker than 
expected based on the strength of its other emission lines, 
notably C IV . These traits have jointly been interpreted as 
evidence for neutral gas in the vicinity of ULAS J1120-I-0641, 
with the weak Lyman-a line being due to damping wing ab¬ 
sorption extending to Arest > 1216 A. In a scenario where 
the neutral gas is intergalactic, models of the quasar’s prox¬ 
imity zone suggest that the surrounding IGM may be >10 
per cent neutral (Bolton et al. 20111. The lack of interven¬ 
ing metal absorption lines further suggests that any neutral 
gas would have to be highly metal poor, particularly if it is 
confined to a discrete absorber (Simcoe et al. 20121. 


The case for neutral gas in the vicinity of ULAS 
J1120-I-0641, whether localized or intergalactic, depends 
strongly on the reality of the damping wing absorption. 
A truncated proximity zone, while unusual, could be due 
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Figure 1. C IV Equivalent Width-Blueshift Anomaly diagram 
showing the C IV emission properties of 5207 SDSS quasars in 
black and ULAS J1120+0641 in blue. This includes all DR7 
quasars in the range 2.4<z<4 for which spectra with STN>5 were 
available. Error bar in the bottom right corner is representative 
for points located in the middle of the distribution; scatter in 
extreme values is due to objects whose C IV line is not well mod¬ 
elled by the spline fit, such as BAL quasars. Thick lines: initial 
cuts from which the matching (blue) and offset (red) selections 
were extracted (see text). The location of the composite quasar 
spectrum from [Mortlock et al.| ( |^011| | is shown in light blue. 
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to an optically thick (A^hi ^ 10^^ cm“^) absorber in the 
vicinity of the qnasar, but a damping wing requires a large 
column density (A^hi sJ cm“^) of neutral gas. A chal¬ 
lenge in analysing such high-redshift quasars, however, is 
that their intrinsic spectra, including their Lyman-a emis¬ 
sion lines, exhibit significant variation between objects. This 
combined with the fact that the blue side of Lyman-a is of¬ 
ten strongly affected by Lyman-a forest absorption (gener¬ 
ally arising from ionized gas) means that the intrinsic shape 
of Lyman-a is largely unknown a priori. This makes it dif¬ 
ficult to estimate the extent to which the emission line flrrx 
has been absorbed, and in the case of ULAS J1120-I-0641, 
whether damping-wing absorption may be present. 

One approach to resolving this issue has been to es¬ 
timate the intrinsic quasar spectrum near Lyman-a using 
either Principal Component Analysis (PCA, [Sukuzi et al.[ 


2005 Paris et al. 20111 or composites made from lower- 


redshift objects (e.g., Vanden Berk et al.| 2001 Cool et 
El |2006[ ). The general approach is to find lower-redshift 
quasars whose spectra match the object of interest redward 
of Lyman-a. The spectra of the lower-redshift objects, which 
are less affected by Lyman-a forest absorption, can then 
be used to predict the unabsorbed continuum around the 
Lyman-a emission line and over the forest. This was the ap¬ 
proach taken by recent studies of ULAS J1120-I-0641 (|Mort-| 


lock et al. 2011 Simcoe et al. 20121. In these cases the com¬ 


posite was constructed primarily to match the C IV emission 
line, whose properties are known to correlate with those of 
Lyman-a (Richards et al. 20111. 

The fidelity of a composite match will naturally depend 
on the availability of similar objects at lower redshifts. How¬ 
ever, the C IV line of quasars exhibits a source systematic ve¬ 
locity shift (Richards et al. 2002 Shang et al. 20071 which 


Selection EWciv/UW^^yAvciv — 

Matching 0.8 —> 1.2 A -300 —^ 300 km 

Offset 0.8 1.2 A 300 900 km s"! 

Table 1. Selection criteria for our initial C IV emission line cuts. 
The matching selection is chosen to match the C IV emission of 
ULAS J1120-I-0641 as well as possible. The offset selection is cho¬ 
sen to mismatch ULAS J1120-|-0641’s C IV emission line blueshift 
by 600 km s“^ or about 3 A. 


is usually attributed to the presence of strong winds and 


jets associated with the quasar (Proga, Stone & Kallman 


20001. In addition the Equivalent Width (EW) of the C IV 


line is known to correlate with Equivalent Width of Lyman- 
a. Mortlock et al. (20111 match ULAS J1120-I-0641 by se¬ 


lecting objects based on their C IV EW and blueshift, but 
acknowledge that this is complicated by the fact that ULAS 
J1120-1-0641 has a large C IV blueshift and is matched by 
relatively few lower-z objects. The sparsity of suitable ob¬ 
jects is aggravated by the need to restrict the selection to 
objects within the narrow redshift range 2.3 < z <2.6 cho¬ 
sen to minimize the impact of Lyman-a forest absorption. 
A compromise thus had to be made in order to obtain a 
large enough sample of objects, resulting in a small misalign¬ 
ment in C IV blueshift between the composite and ULAS 
J1120-|-0641’s spectrum. In this paper we examine whether 
the difference in blueshift, although small, might be hav¬ 
ing an unanticipated effect on the predicted strength of the 
Lyman-a emission line. 


2 METHODS 


We aim to predict the strength of the Lyman-a emission line 
without making use of a composite spectrum. Our approach 
is to select two samples of lower redshift quasars solely on 
their C IV emission properties, regardless of how rare they 
may be. One sample is designed to contain objects simi¬ 
lar to ULAS J1120-I-0641 in C IV properties, except with 
a small blueshift mismatch; the objects in the second sam¬ 
ple are chosen to display no such offset. We then measure 
the Lyman-a flux strength for all objects, which allows us 
to establish the rarity of quasars with Lyman-a emission in¬ 
trinsically as weak as in ULAS J1120-1-0641 in both samples. 

Our sample of comparison spectra is drawn from the 
Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) 
quasar catalog (York et al.l 20001. The DR7 catalog con¬ 


tains 10871 objects in the range 2.4 < z < 4, where both 
the Lyman-a and C IV emission lines are in the observed 
wavelength range. DR7 was chosen because of the availabil¬ 
ity of high quality redshifts measured by [Hewett fc Wild] 
(20101. In this work we use only objects where the systemic 


redshift was determined primarily from the C ill] emission 
line complex. These redshifts were used to shift the spectra 
into the rest frame, and for assessing the relative blueshift 
of the C IV line. 


Similarly to Mortlock et al. (20111, we selected sub¬ 


samples of quasars based on their C IV rest-frame equiv¬ 
alent width, Wciv, and velocity shift relative to systemic. 
Awe IV. We do not attempt to provide an absolute value 
of ULAS J1120-|-0641’s C IV line blueshift, as doing so re¬ 
quires more careful treatment of the underlying continuum 
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Figure 2. Black: C IV emission line of ULAS J1120+0641 from 
GNIRS spectrum. Red: Composite quasar spectrum from [Mort-| 
[lock et al.| 1 |201T| |. Blue: Spline fit to ULAS J1120+0641’s C IV 
emission line. The distance between the peaks of the fits is ~3 A 
in the rest frame of the quasar, corresponding to a velocity offset 
of ~600 kms“^. 


as well as nearby emission lines, in 

larticular Fell. For an 

up-to-date value, see De Rosa et al. 

(20141. Since the ulti- 

mate goal is to identify objects with properties similar to 


ULAS J1120+0641, we used the ULAS J1120+0641 C IV 
emission line itself as a basis for the classification. Using 
the FORS2+GNIRS spectrum presented by |Mortlock et ah] 


(20111, we fit the ULAS J1120+0641 C IV line over the 
wavelength range 1435 < Arest < 1640 A using two com¬ 
ponents, a power law fit over 1435 < Arest < 1480 A and 
1580 < Arest < 1640 A, plus a spline, fit by hand, to the 
emission line flux in excess of the power law shown in Fig. 
2. We then fit each object in the DR7 catalog over the same 
wavelength range using a power law plus an emission line 
template based on the ULAS J1120+0641 spline. The emis¬ 
sion line fit, normalized by the power law continuum, is given 
by 


Ac IV(A') = aA(i\"v°(A), 


( 1 ) 


where Ac\^ is the ULAS J1120+0641 emission line fit, nor¬ 
malized by its power law, and a is a scaling factor. The 
rest-frame wavelengths are related as 


A' = A 1+^^ 




1 + s 


(A - 1540 A) 


A 


( 2 ) 


where Sv is the velocity shift of the C IV line relative 
to ULAS J1120+0641, and s is a stretch factor about 
1540 A, which lines near the peak of the C IV line in ULAS 
J1120+0641. Nominally, objects with C IV lines that are well 
matched to ULAS J1120+0641 will have a ~ 1, ~ 0, and 

s ~ 0. Given these parameters, VFciv for each object as a 
fraction of that of ULAS J1120+0641 can then be recovered 
by Wciv/Whf^ = a(l + s). 

Fig. 1 shows the resulting G IV emission parameter 
space for 7025 SDSS quasars over 2.4 < 2 : < 4.0 for 


which spectra with Hewett & Wild (20101 C ill] redshifts 
were available. As noted by Mortlock et al.| ( |2011 1, the 
ULAS J1120+0641 G IV emission line appears to be highly 
blueshifted compared to the general population of quasars. 


This makes it fall at the edge of the usual distribution of 
CIV parameters. 

We proceeded to select one sample of objects located 
close to ULAS J1120+0641 in C IV parameter space, and an¬ 
other sample with the C IV blueshift offset from the nominal 
value by ~3 A (~600 kms“^). This offset is chosen to mimic 
the apparent offset in C IV between the ULAS J1120+0641 
spectrum and the composite fit of Mortlock et al. (20111, as 
shown in Fig. 2. The cuts we used are shown in Table 1. 

A cut of SNR >5 per pixel, measured by making use 
of the associated SDSS error arrays over the wavelength 
range 1450 - 1500 A, was also applied in the interest of mea¬ 
surement accuracy. For comparison the GNIRS spectrum 
of ULAS J1120+0641 has SNRjii 2 o ~14 per velocity bin 
over the same range. This initial cut yielded unrefined sam¬ 
ples consisting of 204 C IV-matched quasars and 709 offset 
quasars. Due to the large uncertainties in the values of EW 
and blueshift for the objects in our sample, we expected the 
initial cuts to be inefficient at finding all of, and only, ex¬ 
cellent matches for ULAS J1120+0641’s G IV emission. For 
this reason the initial selection boxes are large and a second, 
manual cut was necessary to refine the selection. 

Next, each object’s G IV emission line was overlayed 
with ULAS J1120+0641’s G IV emission line and similarity 
was evaluated visually. Objects which displayed too much 
or too little line asymmetry, as well as those whose C IV 
emission peak was visually offset from ULAS J1120+0641’s 
corresponding value, were rejected and excluded from the 
sample. This yielded a refined sample of 111 objects which 
closely matched ULAS J1120+0641’s G IV emission. The 
same process was repeated for the C IV-offset sample, yield¬ 
ing a refined sample of 216 objects with a C IV blueshift off¬ 
set. The mean redshifts for the matched and offset selections 
were z=2.877 and z=2.912 respectively while the median val¬ 
ues were z=2.898 and z=2.937 respectively. A two-sided KS 
test revealed no evidence that the two samples had been 
drawn from different distributions (p>0.10). 

Measuring the intrinsic Lyman-a flux is particularly dif¬ 
ficult for our samples since they include quasars with red¬ 
shifts as high as z=4, for which Lyman-a forest absorption is 
a greater issue. This is further complicated by the complex 
shape of the Lyman-a emission line over this range, which 
includes contributions from asymmetric Lyman-a emission 
as well as NV. For these reasons the mean unabsorbed con¬ 
tinuum in each object was estimated over the range 1216- 
1220 A, which was chosen to maximize our sensitivity to 
the Lyman-a emission strength, while minimizing the im¬ 
pact of Lyman-a forest absorption. Even in this wavelength 
range, however, Lyman-a and metal absorption lines can 
potentially obscure the intrinsic continuum, particularly at 
higher redshifts. We therefore inspect each object visually 
and estimate the continuum by eye, accounting for absorp¬ 
tion lines to the extend possible. This will naturally be an 
imperfect process; however, the errors in the continuum es¬ 
timate arising from absorption lines are expected to impact 
both samples equally. To mitigate any bias caused by manual 
measurement, we measured the two samples together with 
the objects in random order. As a further check, we calculate 
the mean spectrum of the two samples without attempting 
to correct for absorption lines (see below). 

Emission line fluxes for each object were calculated as 
a fraction of the underlying power-law continuum. The flux 
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normalised CIV peak flux (arb. units) 


Figure 3. Comparison of C IV peak flux and Lyman-a red wing 
flux for 216 quasars matching ULAS J1120+0641 in C IV EW 
and offset from ULAS J1120+0641’s correct CIV Blueshift value 
by ~ 600 kms“^, drawn from the red box in Fig. 1. ULAS 
J1120+0641’s location is shown by a thick blue star. Red star 
indicates the location of the [Mortlock et al!] l |2011[ | composite 
quasar spectrum when measured in the same way. The Lyman-a 
red wing flux is measured between 1216A- 1220A. Among this 
population ULAS J1120+0641 is a 97 per cent outlier in Lyman- 
a red wing flux. Representative error bar shown in the bottom 
right corner. 


of each object was normalised by measuring the mean flux 
over the wavelength range 1450 - 1500 A and the continuum 
emission was modeled by a power-law fitted iteratively over 
narrow, emission line free regions of the spectrum. These 
fitting windows spanned the rest wavelength ranges 1320 to 
1350 A, 1435 to 1480 A and 1600 to 1700 A. The range 2000 
- 2050 A was nsed in addition for z<3 objects for which 
those wavelengths were present in the rest frame. The peak 
flux over C IV was measured by the same method for all 
objects as well as ULAS J1120+0641 to check that the C IV 
emission was well matched by our selections and to ensure 
that direct comparison with the flux value over the Lyman-a 
red wing was meaningful. 


3 RESULTS AND DISCUSSION 


Among the objects which matched the C IV emission in 
shape and EW but which were offset in blueshift by ~ 600 
km s“^, only three per cent of objects showed Lyman-a emis¬ 
sion weaker than ULAS J1120+0641 over the same range 


(Fig. 3). This is consistent with the assessment by Mort¬ 


lock et al. (20111 that ULAS J1120+0641 is a ~2cr out¬ 


lier from the composite they used. In contrast, among the 
selection which matched ULAS J1120+0641 in C IV emis- 


Figure 4. Comparison of C IV peak flux and Lyman-a red wing 
flux for 111 quasars matching ULAS J1120+0641’s C IV emission 
in both EW and Blueshift, drawn from the blue box in Fig. 1. 
ULAS J1120+0641’s location is shown by a thick blue star. The 
objects highlighted in green are shown in Fig. 6. The Lyman-a 
red wing flux is measured between 1216A- 1220A. Among this 
population, 31 per cent of objects have Lyman-a lines that are 
weaker than ULAS J1120+0641. Representative error bar shown 
in the bottom right corner. 


sion, we find that 31 per cent have Lyman-a emission over 
the red wing weaker in magnitude than ULAS J1120+0641 
(Fig. 4). In addition our C IV measurements are shown to 
demonstrate no residual correlation in these small samples 


(Rmatched = 0.20, Roffaet = 0.07), demonstrating that the 
C IV emission line of ULAS J1120+0641 is typical among 
these samples. We therefore find no compelling evidence that 
the Lyman-a flux of ULAS J1120+0641 is anomalously low. 

The averages and 68 (95) per cent bounds in Lyman-a 
flux values for the matched and offset sample are Flc = 
2-64lgJ| (toito) and Flc = 2.88lJ5;J^ (ias”) respectively, 
corresponding to variations of ~ ±10 (~ ± 20) per cent 
in both distributions. These values are similar to the cor¬ 
responding value of 13 per cent quoted by [Mortlock et al.| 
(20111 for the spread in Lyman-a flux once C IV blueshift 
is constrained; however the distributions are mildly non- 
gaussian since they have a tail and those numbers are given 
as a consistency check only. In addition these values of the 
spread are given at a particular wavelength. 

Small sample size, low SNR and scatter mean that our 
sample might not be suitable for making an appropriate 
composite fit to ULAS J1120+0641, and such an analysis 
is beyond the scope of this work. Fig. 5 presents the aver¬ 
age spectrum for both our selections; those were made by 
interpolating the spectra onto a common wavelength array 
and taking the mean value pixel-by-pixel. The spectra were 
power-law corrected and normalised but absorption was not 
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Figure 5. Average spectra of the matched selection of quasars (black) and the offset selection (red). Note how a small mismatch in 
the C IV emission line leads to large variation of the Lyman-o peak flux, while the low ionisation lines are not affected. The spectra 
were normalized by dividing by a fitted power-law; absorption was not taken into account and as a consequence these averages do not 
represent the underlying continuum emission. Uncorrected absorption tends to both lower and smooth the peaks hence the difference in 
mean spectra is less pronounced than in the ensemble distributions (see Fig. 3 and 4). 


addressed. The difference in flux between the peaks of the 
average spectra in Fig. 5 is comparable to the difference 
between the averages of the distributions in Fig. 3 and 4; 
however the effect of absorption is to weaken the emission 
feature in the average spectra. This confirms the fact that 
a small mismatch in the C IV emission line can have a sig¬ 
nificant effect on Lyman-a emission. Although this effect is 
not large, it is sufficient to significantly impact the degree 
to which ULAS J1120-1-0641 is an outlier. The amount by 
which the average spectra differ in C IV blueshift is similar to 
the mismatch between ULAS J1120-1-0641 and its composite 
fit as published in |Mortlock et aL] ( |2011[ ), condoning our use 
of the value of a ~ 600 kms“^ C IV blueshift mismatch. 


Example objects matching ULAS J1120+0641’s spec¬ 
trum extremely well over its entire continuum and down 
to the onset of the Lyman-a forest at 1216A are shown 
in Fig. 6 along with their respective redshifts. Using the 
strength of flux right of 1216A in a manner similar to ULAS 
J1120-I-0641, one would thus falsely infer the presence of a 
damping wing in those objects. A large variety of Lyman-a 
peak shapes exists among the objects which match ULAS 
J1120-|-0641’s C IV line; in particular the N V emission line 
at 1240A overlaps with Lyman-a in a non-trivial way es¬ 
pecially when the latter is weak. We make no attempt to 
constrain the N V emission. The majority of objects in our 
sample display a visible Si II emission line at 1260A, suggest¬ 
ing that the high-ionization lines C IV , N V and Lyman-o? 
are blueshifted together while low-ionization lines Si II , Ol, 
CII are not. Our findings are in agreement with Richards et 
al. (2011) who identify a link between a large NV/Lyman-a 
ratio in quasars with extreme C IV blueshifts as well as a 
possible C IV blueshift/Lyman-a blueshift correlation. This 
effect also helps to resolve the apparent mismatch of N V in 
the Mortlock et al.l (20111 composite. 


4 CAVEATS 


We have argued that the scatter in Lyman-a flux for a 
given C IV EW, combined with a small decrease in Lyman-a 
flux with increasing C IV blueshift, is sufficient to bring the 
red wing of the Lyman-a line of ULAS J1120-1-0641 within 
the distribution seen in lower-redshift objects. Reconstruct¬ 
ing the unabsorbed continuum of a quasar is known to be 


a challenging problem, however (e.g., Paris et al. 2011 and 
references therein), and it is important to bear in mind cer¬ 
tain caveats. For example, accurate measurements of C IV 
blueshifts require reliable estimates of the systemic quasar 
redshifts. For the comparison SDSS spectra we have used 
redshifts based on C ill] from |Hewett &; Wild| | [20To| ) , assum¬ 
ing that these are independent from the C IV and Lyman-a 
properties. The C ill] feature is a blend of Al ill, Si ill, and 
C III] emission, however (e.g., [Vanden Berk et al.[ |200ip , 
and so the measured redshift will depend on the relative 
strengths and velocity shifts of these components. If the 
properties of this emission complex correlate with the prop¬ 
erties of the C IV and/or Lyman-a lines, then this may intro¬ 
duce systematic biases into the relationship we infer between 
C IV and Lyman-a. 

As noted earlier, the ULAS J1120-1-0641 redshift is mea¬ 
sured from [C ll] 158 pm emission from the ISM of the 
host galaxy ( [Venemans et ah] [2012 1. There is some indi¬ 
cation, however, that C ill] in ULAS J1120-I-0641 is slightly 
blueshifted with respect to [C ll] (see Fig. 1 of Mortlock 
et al. 20111. If this is generally true for quasars with large 


C IV blueshifts, then we should potentially adopt the C ill] 
redshift for ULAS J1120-1-0641 when comparing it to objects 
with redshifts derived from C ill]. This would give a smaller 
relative blueshift for C IV in ULAS J1120-I-0641, meaning 
that our comparison objects should be drawn from a re¬ 
gion closer to the red box in Fig. 1. By the arguments pre¬ 
sented here, this would imply a somewhat higher expected 
intrinsic Lyman-a flux over 1216-1220 A. On the other hand. 
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Figure 6. SDSS spectra of objects at lower redshifts which match J1120+0641’s spectrum in CIV as well all the way to the onset of 
the Lyman-a forest at 1216A, indicated by the blue dashed line. The GNIRS spectrum of ULAS J1120+0641 , binned in bins of 3, is 
overlayed in red. All spectra are normalized by a fitted power-law. SDSS spectra are binned in bins of 5. The C IV parameters of these 
objects are shown as green dots on Fig. 4. 
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the observed wavelength at which discrete Lyman-a forest 
absorption begins in ULAS J1120+0641, i.e., the shortest 
wavelength at which a smooth damping could be detected, 
would correspond to a longer rest-frame wavelength, and 
hence a weaker portion of the Lyman-a emission line pro¬ 
file. It is therefore unclear whether adopting a lower redshift 
for ULAS J1120-I-0641 would strengthen the case for damp¬ 
ing wing absorption. 

A robust reconstruction of the intrinsic Lyman-a flrrx 
may ultimately require a more sophisticated analysis than 
the one presented here. We have shown that there is some 
uncertainty in the Lyman-a prohle related to the C IV prop¬ 
erties, and hence motivation to re-examine the claim of a 
damping wing. A conclusive analysis, however, may need 
to take into account the detailed properties of all available 
emission lines, and be verified through tests on lower-redshift 
quasars with known Lyman-a profiles. 

5 CONCLUSIONS 


Composites quasar spectra have been used in the past to 
study the IGM around ULAS J1120-I-0641 and hypothesize 
a large fraction of neutral gas or very metal-poor gas. Those 
claims were based on ULAS J1120+0641’s seemingly highly- 
absorbed Lyman-a emission (Mortlock et al. (20111 using 
the quasar composite technique from Hewett & Wild (20101; 


Simcoe et al. (20121; Bolton et al. (20111). However those 


authors acknowledge the difficulty of producing a composite 
ht to ULAS J1120-I-0641, due primarily to the rarity of suit¬ 
able objects. As a result, composite hts have mismatched its 
C IV emission line by ~ 600 kms“^ or more in blueshift. 


Mortlock et al. (20111 suggested that this mismatch might 


have an impact on the existence of a Lyman-a damping 
wing. 

We have tested the effect of this small C IV mismatch by 
selecting a sample of 111 SDSS DR7 quasars at 2.4 < 2 < 4 
that more closely match ULAS J1120-|-0641’s C IV emission, 
and a second sample of 216 quasars whose C IV emission 
lines match ULAS J1120-|-0641’s in EW but are offset blue- 
ward by ~600 kms”'^ (~3A in the rest frame). We find that 
among a population of quasars which match its C IV emis¬ 
sion in equivalent width and shape, but offset by ~ 600 
kms“^ from the correct C IV line systematic blueshift, 
ULAS J1120-1-0641 appears anomalous to a confidence of 
97 per cent, in agreement with Mortlock et al. (20111 who 
find ULAS J1120+0641 to be a 2a outlier from their com¬ 
posite fit. When compared to a population of lower-redshift 
quasars which match its C IV line emission, however, ULAS 
J1120+0641’s Lyman-a flux does not appear to be anoma¬ 
lously weak. 

In light of our results we suggest that the shape of the 
Lyman-a line in ULAS J1120+0641 is not obviously due 
to a damping wing. We note that the case for neutral gas 
in the vicinity of ULAS J1120+0641 depended on both the 
damping wing absorption and the shortness of the prox¬ 
imity zone ( [Bolton et al.| |201ip . The latter could be due 
to an optically thick absorber near the quasar redshift, a 
point we will address in a future paper. Independent of the 
ULAS J1120+0641 constraints, however, there is evidence 
for a significantly neutral IGM at redshifts z>7 from the 
evolution of the Lyman-a emission fraction in galaxies (e.g., 
Caruana et al.| |2012| [Treu et ah] |2013| [Tilvi et al.[ |2014[ 
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quasars would signihcantly strengthen this claim. Larger 
samples of quasars at z>7, along with a more robust method 
to reconstruct their intrinsic Lyman-a flux, are therefore of 
great interest for future reionisation studies. 
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